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ABSTRACT 

We use a large sample of isolated dark matter halo pairs drawn from cosmological 
N-body simulations to identify candidate systems whose kinematics match that of the 
Local Group of Galaxies (LG). We find, in agreement with the “timing argument” and 
earlier work, that the separation and approach velocity of the Milky Way (MW) and 
Andromeda (M31) galaxies favour a total mass for the pair of ~ 5 x 10^^ M 0 . A mass 
this large, however, is difficult to reconcile with the small relative tangential velocity 
of the pair, as well as with the small deceleration from the Hubble flow observed for 
the most distant LG members. Halo pairs that match these three criteria have average 
masses a factor of ^ 2 times smaller than suggested by the timing argument, but with 
large dispersion, spanning more than a decade in mass. Guided by these results, we 
have selected 12 halo pairs with total mass in the range 1.6-3 .6 x 10^^ Mq for the 
APOSTLE project (A Project Of Simulations of The Local Environment), a suite of 
resimulations at various numerical resolution levels (reaching up to 10 ^ M 0 per gas 
particle) that use the hydrodynamical code and subgrid physics developed for the EA¬ 
GLE project. These simulations reproduce, by construction, the main kinematics of the 
MW-M31 pair, and produce satellite populations whose overall number, luminosities, 
and kinematics are in good agreement with observations of the MW and M31 com¬ 
panions. These diagnostics are sensitive to the total mass assumed for the MW-M31 
pair; indeed, the LG satellite population would be quite difficult to reproduce for pair 
masses as high as indicated by the timing argument. The APOSTLE candidate systems 
thus provide an excellent testbed to confront directly many of the predictions of the 
ACDM cosmology with observations of our local Universe. 

Key words: Cosmology, Numerical simulations, Dark matter, Local Group, dwarf 
galaxies 


1 INTRODUCTION 


for studies of the formation and evolution of galaxies. Their 
close vicinity implies that LG galaxies are readily resolved 
into individual stars, enabling detailed exploration of the 
star formation, enrichment history, structure, dark matter 
content, and kinematics of systems spanning a wide range 
of masses and morphologies, from the two giant spirals that 
dominate the Local Group gravitationally, to the faintest 
galaxies known. 


The Local Group of galaxies (LG), which denotes the as¬ 
sociation of the Milky Way (MW) and Andromeda (M31), 
their satellites, and galaxies in the surrounding volume out 
to a distance of ~ 3 Mpc, provides a unique environment 


* Email: azadehf@uvic.ca 


2 Fattahi et al. 


This level of detail comes at a price, however. The Local 
Group volume is too small to be cosmologically representa¬ 
tive, and the properties of its galaxy members may very 
well have been biased by the peculiar evolution that led to 
its particular present-day conhguration, in which the MW 
and M31, a pair of luminous spirals ~ 800 kpc apart, are ap¬ 
proaching each other with a radial velocity of 120 km s“^. 
This galaxy pair is surrounded by nearly one hundred galax¬ 
ies brighter than Mv ~ —8, about half of which cluster 
tight ly around our Galaxy and M31 (see, e.g.. lMcGonnachI3 
I201I for a recent review). 

The Local Group is also a relatively isolated environ¬ 
ment whose internal dynamics are dictated largely by the 
MW-M31 pair. Indeed, outside the satellite systems of MW 
and M31, there are no galaxies brighter than Mb = —18 (the 
luminosity of the Large Magellanic Gloud, hereafter LMC 
for short) within 3 Mpc from the MW. The nearest galaxies 
comparable in brightness to the MW or M31 are just beyond 
3.5 Mpc away (NGG 5128 is at 3.6 Mpc; M81 and NGC 253 
are located 3.7 Mpc from the MW). 

Understanding the biases that this particular envi¬ 
ronment may induce on the evolution of LG members 
is best accomplished through detailed numerical simula¬ 
tions that take these constraints directly into account. 
This has been recognized in a number of recent studies, 
which have followed small volumes tailore d to resemble, in 
broad terms, the Local Group (see , e.g., iGottloeber et all 
I2OI0I : iGarrison-Kimmel et al.ll2014h . This typically means 
selecting ~ 3 Mpc-radius regions where the mass bud¬ 
get is dominated by a pair of virialized halos separated 
by the observed MW-M31 distance and whose masses are 
chosen to match various additional constraints (see, e.g.. 


Forero-Romero et al.ll2013ll. 

The mass constraints may include estimates 
individual viria|3 masses of both MW and M31, 
cally based on the kinematics of tracers such as 
lite galaxies, halo stars, or tidal debris (see. 

jf the 
typi- 
satel- 

e.g., 

Battaglia et al.l l2005l; Sales et al.l 20071: Smith et al.j 

20071: 

Xue et al.l l2008l: IWatkins et al.l I2OI0I: iDeason et al.l 

20121: 

Bovlan-Kolchin et al.l l2013l: IPifll et al.l I20l4 iBarber et al.l 


2014 for some recent studies). However, these estimates are 


usually accurate only for the mass enclosed within the re¬ 
gion that contains each of the tracers, so that virial mass 
estimates are subject to non-negligible, and potentially un¬ 
certain, extrapolation. 

Alternatively, the MW and M31 stellar masses may be 
combined with “abun dance matching” techniques to derive 
virial masses (see, e.g ., IGuo et ^l2ninl : iBehroozi et al.l[2013l : 
iKravtsov et al.l[20l4 and references therein). In this proce¬ 
dure, galaxies of given stellar mass are assigned the virial 
mass of dark matter halos of matching number density, com¬ 
puted in a given cosmological model. Shortcomings of this 
method include its reliance on the relative ranking of halo 
and galaxy mass in a particular cosmology, as well as the 
assumption that the MW and M31 are average tracers of 
the halo mass-galaxy mass relation. 


^ We define the virial mass, M 200 , as that enclosed by a sphere 
of mean density 200 times the critical density of the Universe, 
Pcrit = /8nG. Virial quantities are defined at that radius, 

and are identified by a “200” subscript. 


A further alternative is to use the kinematics of LG 
members to estimate virial masses. On e example is the “nu - 
merical action” method developed by IPeebles et al.l (|20nih 
to reconstruct the peculiar velocities of nearby galaxies 
which, when applied to th e LG, predicts a fai rly large cir¬ 
cular velocity for the MW JPeebles et al.ll2nilh . A simpler, 
but nonetheles s useful, example is p rovided by the “tim¬ 
ing argument” llKahn fc Woltieilll959l l. where the MW-M31 
system is approximated as a pair of isolated point masses 
that expand radially away after the Big Bang but decel¬ 
erate under their own gravity until they turn around and 
start approaching. Assuming that the age of the Universe is 
known, that the orbit is strictly radial, and that the pair is 
on first approach, this argument leads to a robu st and unbi¬ 
ased estimate of the total mass of the system fc Whitel 
I2OO8II . Difficulties with this approach include the fact that 
the tangential veloci ty of the pair is neglected (see, e.g., 
iGonzalez et al.ll2014h . together with uncertainties relating 
the total mass of the point-mass pair to the virial masses of 
the individual systems. 

Finally, one may use the kinematics of the outer LG 
members to estimate the total mass of the MW-M31 pair, 
since the higher the mass, the more strongly LG members 
should have been decelerated from the Hubble flow. This 
procedure is appealing because of its simplicity but suffers 
from the uncertain effects of nearby massive structures, as 
well as from difficulties in accounting for the directional de¬ 
pendence of the deceleration and, possibly, for the gravita¬ 
tional tor que/pull of even more distant large-scale structure 
(see, e.g.. |Penarrubia et al.ll20l4lSorce et al.ll2of3 . for some 
recent work on the topic). 

A review of the literature cited above shows that these 
methods produce a range of estimates (spanning a factor 
of 2 to 3) of the individual masses of the M W and M31 
and/ or the total mass of the Local Group f see I Wang et al.l 
l2015l. for a recent compilation). This severely conditions 
the selection of candidate Local Group environments that 
may be targeted for resimulation, and is a basic source of 
uncertainty in the predictive ability of such simulations. 
Indeed, varying the mass of the MW halo by a factor 
of 3, for example, would likely lead to variations of the 


such systems 

IBov 

an-Kolchin et al.ll201ll:IWang et al.ll2012l: 

ICautun et al.ll2014 

), limiting the insight that may be gained 


from direct quantitative comparison between simulations 
and observations of the Local Group. 

With these caveats in mind, this paper describes the se¬ 
lection procedure, from a simulation of a large cosmological 
volume, of 12 viable Local Group environment candidates 
for resimulation. These 12 candidate systems form the basis 
of the EAGLE-APOSTLE project, a suite of high-resolution 
cosmological hydrodynamical resimulations of the LG envi¬ 
ronment in the ACDM cosmogony. The goal of this paper 
is to motivate the particular choices made for this selec¬ 
tion whilst critically reviewing the constraints on the total 
mass of the Local Group placed by the kinematics of LG 
members. Preliminary results from the project (which we 
shall hereafter refer to as APOSTLE, a shorthand for “A 
Project Of Simulations of The Local Environment”), which 
uses the same code developed for t he EAGLE simulations 
l|Schave et al.ll2015l:lGrain et al.ll2015lL have already been re ¬ 
ported in ISawala et al.l ll20l4l2015h and lOman et ah limi). 
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Table 1. The parameters of the cosmological simulations used in this paper. 


Simulation 

Cosmology 


Ha 

Hf, 

h 

0-8 

ris 

Cube side 
(Mpc) 

Particle number rrip dm 

(Mo) 

MS-I 

WMAP-l 

0.25 

0.75 

0.045 

0.73 

0.9 

1 

685 

2160® 

1.2 X 10® 

MS-II 

WMAP-1 

0.25 

0.75 

0.045 

0.73 

0.9 

1 

137 

2160® 

9.4 X 10® 

EAGLE (L100N1504) 

Planck 

0.307 

0.693 

0.04825 

0.6777 

0.8288 

0.9611 

100 

15043 

9.7 X 10® 

DOVE 

WMAP-7 

0.272 

0.728 

0.0455 

0.704 

0.81 

0.967 

100 

1620® 

8.8 X 10® 


We beRin by using the Millenniu m _ Simulations 

(ISpringel et al.l l2005l : iBovlan-Kolchin et al.l l2009l l to select 
relatively isolated halo pairs separated by roughly the dis¬ 
tance between the MW and M31 and derive the distribution 
of total masses of pairs that reproduce, respectively, the rel¬ 
ative radial velocity of the MW-M31 pair, or its tangential 
velocity, or the Hubble flow deceleration of distant LG mem¬ 
bers. 

Given the disparate preferred masses implied by each of 
these criteria when applied individually, we decided to select 
pairs within a narrow range of total mass that match loosely 
the LG kinematics rather than pairs that match strictly the 
kinematic criteria but that span the (very wide) allowed 
range of masses. This choice allows us to explore the “cosmic 
variance” of our results given our choice of LG mass, whilst 
guiding how such results might be scaled to other possible 
choices. We end by assessing the viability of our candidate 
selection by comparing their satellite systems with those of 
the MW and M31 galaxies. 

The plan for this paper is as follows. We begin by assess¬ 
ing in Sec. [2] the constraints on the LG mass placed by the 
kinematics of the MW-M31 pair and other LG members. We 
describe next, in Sec. O the choice of APOSTLE candidates 
and the numerical resimulation procedure. Sec. |4] analyzes 
the properties of the satellite systems of the main galaxies 
of the LG resimulations and compares them with observed 
LG properties. We end with a brief summary of our main 
conclusions in Sec. [5] 


2 THE MASS OF THE LOCAL GROUP 
2.1 Observational Data 


We use below the positions, Galactocentric distances, line- 
of-sight-velocities and V-band magnitudes (converted to 
stellar masses assuming a mass-to-light ratio of unity in so¬ 
lar units) of Local Grou p mem bers as given in the com¬ 
pilation of iMcConnachij (l2012l i. We also use the relative 
tangential velocit y of the M31-MW pair der ived from M31’s 
proper motion by van der Marel et al.l ll2012l ~). When needed, 
we assume an LSR velocity of 220 km s“^ at a distance of 
8.5 kpc from the Galactic center and that the Sun’s peculiar 
motion relative to t he LSR is Up) = 11.1, Uq = 12.24 and 
Wq = 7.25 km s“^ dSchonrich et al.ll2010l '). to refer veloci¬ 
ties and coordinates to a Galactocentric reference frame. 

According to these data, the MW-M31 pair is 787 ± 
25 kpc apart, and is approaching with a relative radial ve¬ 
locity of 123 ± 4 km s“^. In comparison, its tangential ve¬ 
locity is quite low: only 7 km s“^ with Icr confidence region 
Si 22 km s“^ . We shall assume hereafter that these values 


are comparable to the relative velocity of the centers of mass 
of each member of the pairs selected from cosmological sim¬ 
ulations. In other words, we shall ignore the possibility that 
the observed relative motion of the MW-M31 pair may be 
affected by the gravitational pull of their massive satellites; 
i.e., the Magellanic Clouds (in the case of the MW) and/or 
M33 (in the case of M31). This choice is borne out of simplic¬ 
ity; correcting for the possible displacement caused by these 
massive satellites requires detailed assumptions about their 
orbits an d their masses, whic h are fairly poorly constrained 
(see, e.g.. lG6mez et al.]l2015l '). 

We shall also consider the recession velocity of distant 
LG members, measured in the Galactocentric frame. This 
is also done for simplicity, since velocities in that frame are 
more straightforward to compare wi th velocities measured in 
simu lations. Other work (see, e.g., iGarrison-Kimmel et al.l 
l2014ll has used velocitie s expressed in the Local Group - 
centric frame defined by [Karachentsev fc Makaro"^ (Il996ll . 
This transformation aims to take into account the apex of 
the Galactic motion relative to the nearby galaxies in order 
to minimize the dispersion in the local Hubble flow. This 
correction, however, is sensitive to the volume chosen to 
compute the apex, and difficult to replicate in simulations. 


2.2 Halo pairs from cosmological simulations 

We u se the Mille nnium Simulations. MS-I llSpringel et al.l 
l2005li and MS-II llBovlan-Kolchin et al.|[^009ll . to search for 
halo pairs with kinematic properties similar to the MW and 
M31. The MS-I run evolved 2160® dark matter particles, 
each of 1.2 x 10® M©, in a box 685 Mpc on a side adopt¬ 
ing AGDM cosmological parameters consistent with the 
WMAP-1 measurements. The MS-H run evolved a smaller 
volume (137 Mpc on a side) using the same cosmology and 
number of particles as MS-I. Each MS-H particle has a mass 
of 9.4 X 10® Mq. We list in Table [T] the main cosmological 
and numerical parameters of the cosmological simulations 
used in our analysis. 

At 2 = 0, dark matter halos in bot h simulations were 
identified using a friends-of-friends (FoF, iDavis et al.lll985l 'l 
algorithm run with a linking length equal to 0.2 times 
the mean interparticle separation. Each FoF halo was then 
searched iteratively for self-boun d substructures (sub halos) 
using the SUBFIND algorithm (ISpringel et al.l [20011 1. Our 
search for halo pairs include all pairs of separate FoF halos, 
as well as single FoF halos with a pair of massive subhalos 
satisfying the kinematic and mass conditions we list below. 
The latter is an important part of our search algorithm, since 
many LG candidates are close enough to be subsumed into 
a single FoF halo at 2 = 0. 
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Figure 1. Left-. The relative radial velocity vs distance of all halo pairs selected at highest isolation (“Hilso”) from the Millennium 
Simulation. Colours denote the total mass of the pairs (i.e., the sum of the two virial masses), as indicated by the colour bar. The starred 
symbol indicates the position of the MW-M31 pair in this plane. The dotted line illustrates the evolution of a point-mass pair of total 
mass ~ 5 X 10^^ Mq in this plane (in physical coordinates) which, according to the timing argument, ends at M31’s position. Dashed 
lines indicate the loci, at z = 0, of pairs of given total mass (as labelled) but different initial energies, according to the timing argument. 
The box surrounding the M31 point indicates the range of distances and velocities used to select pairs for further analysis. Right: Same 
as left panel, but for the relative tangential velocity. Dashed curves in this case indicate the mean distance-velocity relation for pairs of 
given total mass, as labelled. 


The list of MS-I and MS-II halos retained for analysis 
include all pairs separated by 400 kpc to 1.2 Mpc whose 
members have virial masses exceeding 10^^ M© each but 
whose combined mass does not exceed 10^® M©. These pairs 
are further required to satisfy a fiducial isolation criterion, 
namely that no other halo more massive than the less mas¬ 
sive member of the pair be found within 2.5!Mpc from the 
center of the pair (we refer to this as “medium isolation” 
or “Medlso”, for short). We have also experimented with 
tighter/looser isolation criteria, enforcing the above crite¬ 
rion within 1 Mpc (“loosely isolated” pairs; or “Lolso”) or 
5!Mpc (“highly isolated” pairs; “Hilso”). Since the nearest 
galaxy with mass comparable to the Milky Way is located 
at ~ 3.5 Mpc, the fiducial isolation approximates best the 
situation of the Local Group; the other two choices allow 
us to assess the sensitivity of our results to this particular 
choice. 


2.3 Radial velocity constraint and the timing 
argument 

The left panel of Fig. [1] shows the radial velocity vs separa¬ 
tion of all pairs in our MS-I samples, selected using our max¬ 
imum isolation criterion. Each point in this panel is coloured 
by the total mass of the pairs, defined as the sum of the virial 
masses of each member. 

The clear correlation seen between radial velocity and 
mass, at given separation, is the main prediction of the “tim¬ 
ing argument” discussed in Sec. [T] Timing argument pre¬ 
dictions are shown by the dashed lines, which indicate the 


expected relation for pairs with total mass as stated in the 
legend. This panel shows clearly that low-mass pairs as dis¬ 
tant as the MW-M31 pair are still, on average, expanding 
away from each other (positive radial velocity), in agreement 
with the timing argument prediction (top dashed curve). In¬ 
deed, for a total mass as low as 2 x 10^^ Mq, the binding 
energy reaches zero at r = 914 kpc, 14 = 43 km s“^, where 
the dashed line ends. 


It is also clear that predominantly massive pairs 
have approach speeds as large as the MW-M31 pair (~ 
— 120 km s“^). We illustrate this with the dotted line, which 
shows the evolution in the r-I4 plane of a point-mass pair of 
total mass 5 x 10^^ Mq, selected to match the MW-M31 pair 
at the present time. The pair reached “turnaround” (i.e., 
null radial velocity) about 5 Gyr ago and has since been ap¬ 
proaching from a distance of 1.1 Mpc (in physical units) to 
reach the point labelled “M31” by z=0, in the left-panel of 

Fig.m 

An interesting corollary of this observation is that 
the present turnaround radius of the Local Group is ex¬ 
pected to be well beyond 1.1 Mpc. Assuming, for guid- 
ance, that the tur naround radius grows roughly like 
ijBertschingerj flQSsI ) . this would imply a turnaround radius 
today of roughly ~ 1.7 Mpc, so that all LG members just 
inside that radius should be on first approach. We shall re¬ 
turn to this point when we discuss the kinematics of outer 
LG members in Sec. [23] below. 
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Figure 2. Radial velocity vs distance for all Local Group members out to a distance of 3 Mpc from the Galactic center. Each panel 
correspond to pairs of different mass, as given in the legends. Symbols repeat in each panel, and c orresp ond to observations in the 
Galactocentric reference frame, taken from the compilations of I McConnachi j ll2012f) and [Tullv et al.l ll2009[) . Solid symbols are used to 
show LG members that lie, in projection, within 30 degrees from the direction to M31. The dotted line in each panel is the timing 
argument curve for M31, as in Fig.^ The coloured lines and shaded regions correspond, in each panel, to the result of binning Medlso 
MS-II halo pairs separated by 600-1000 kpc and least-square fitting the recession velocities of the outer members (between 1.5 and 3 
Mpc from either primary). The solid coloured line shows the median slope and zero-point of the individual fits. The shaded regions show 
the interquartile range in zero-point velocity. Note that the recession speeds of outer members decrease with increasing LG mass. The 
dot-dashed line shows the unperturbed Hubble flow, for reference. 


2.4 Tangential velocity constraints 

The right-hand panel of Fig. ^ shows a similar exercise to 
that described in Sec. da but using the relative tangen¬ 
tial velocity of the pairs. This is compared with that of the 
M W-M31 pair , whic h is measured to be only ~ 7 km s“^ 
by ISohn et al.l (l2012l i and is shown by the starred symbol 
labelled “MSI”. 

This panel shows that, just like the radial velocity, the 
tangential velocity also scales, at a given separation, with 
the total mass of the pair. In general, higher mass pairs 
have higher speeds, as gleaned from the colours of the points 
and by the three dashed lines, which indicate the average 
velocities of pairs with total mass as labelled. 

The average relative tangential velocity of a 5 x 10^^ Mq 
pair separated by ~ 800 kpc is about ~ 100 km s“^, and 
very few of such pairs (only ~ 6%) have velocities as low 
as that of the MW-M31 pair. The orbit of a typical pair 
of such mass is thus quite different from the strictly radial 
orbit envisioned in timing argument estimates. Indeed, the 
low tangential velocity of the MW-M31 pair clearly favours 
a much lower mass for the pair than derived from the timing 
argument Isee lGonzalez et al.ll2014l . for a similar finding). 

The kinematics of the MW-M31 pair is thus peculiar 
compared with that of halo pairs selected from cosmological 
simulations: its radial velocity is best matched with rela¬ 
tively large masses, whereas its tangential velocity suggests 
a much lower mass. We shall return to this issue in Sec. \m 


after considering next the kinematics of the outer LG mem¬ 
bers. 


2.5 Hubble flow deceleration 


The kinematics of the more distant LG members is also sen¬ 
sitive to the total mass of the MW-M31 pair. In particular, 
we expect that the larger the mass, the more the recession 
velocities of those members would be decelerated from the 
Hubble flow. We explore this in Fig. [2] where the symbols 
in each panel s how the Galactocentr ic radial velocity of all 
galaxies in the iMcConnachid (l2012l~l catalogu e and in the 
Extragalactic Distance Database of Tullv et al.l (l2009fl found 
within 3 Mpc of the MW. These data illustrate two inter¬ 
esting points. One is that all galaxies beyond ~ 1.3 Mpc 
from the MW are receding; and the second is that the mean 
velocity of all receding galaxies is only slightly below the 
Hubble flow, Vr = Hor, indicated by the dot-dashed line for 
Ho = 70.4 km s~^ Mpc“^. 

The hrst result is intriguing, given our argument in 
Sec. 12.31 that the turnaround radius of the LG should be 
around ~ 1.7 Mpc at the present time. The second point is 
also interesting, since it suggests that the local Hubble flow 
around the MW (beyond ~ 1.3 Mpc) has been relatively 
undisturbed. The first point suggests that the M31 motion 
is somewhat peculiar relative to that of the rest of the LG 
members; the second that the total mass of the MW-M31 
pair cannot be too large, for otherwise the recession veloci¬ 
ties would have been decelerated more signihcantly. 
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Figure 3. Parameters of the least-squares fits to the recession 
speeds of outer LG members as a function of distance. The solid 
starred symbol labelled “LG” indicates the result of fitting the 
distance vs Galactocentric radial velocity of all LG galaxies be¬ 
tween 1.5 and 3 Mpc from the MW. The Hubble flow is shown by 
the solid circle. Each coloured point corresponds to a halo pair 
selected from the MS-II simulation assuming medium isolation, 
and uses all halos in the 1.5-3 Mpc range resolved in MS-II with 
more than 100 particles (i.e., masses greater than 1 x 10^ Mq). 
Velocities and distances are measured from either primary, and 
coloured according to the total mass of the pair. Note that re¬ 
cession velocities decrease steadily as the total mass of the pair 
increases. The square surrounding the “LG” point indicates the 
error in the slope and zero-point, computed directly from the fit 
to the 33 LG members with distances between 1.5 and 3 Mpc. 
Open starred symbols correspond to the 12 pairs selected for the 
APOSTLE project (see Sec. I3.2ll . 


We examine this in more detail in the four panels of 
Fig.il where the coloured lines show the mean recession 
speed as a function of distance for halos and subhalos sur¬ 
rounding candidate MS- pairs of different mass, binned 
as listed in the legend. Recession speeds are measured from 
each of the primaries for systems in the distance range 1.5- 
3 Mpc and are least-square fitted independently to derive 
a slope and velocity zero-point for each pair. The coloured 
lines in Fig. i]are drawn using the median slope and zero- 
point for all pairs in each mass bin. The shaded areas indi¬ 
cate the interquartile range in the zero-point velocity. This 
figure shows clearly that the more massive the pair, the 
lower, on average, the recession velocities of surrounding sys¬ 
tems. 

This is also clear from Fig.[3l where we plot the individ¬ 
ual values of the slopes and zero-point velocities (measured 
at r = 2.5 Mpc), compared with the values obtained from a 
least-squares fit to the Local Group data. The zero-point of 


^ We use here only MS-II pairs; the numerical resolution of the 
MS-I simulation is too coarse to identify enough halos and subha¬ 
los around the selected pairs to accurately measure the velocities 
of outer LG members. 


the fit is most sensitive to the total mass of the pair; as a 
result, the relatively large recession velocity of the outer LG 
members clearly favours a low total mass for the MW-M31 
pair. 


2.6 Mass distributions 

The histograms in Fig. 2] summarize the results of the pre¬ 
vious three subsections. Each of the top three panels shows 
the distribution of the tatal masses of pairs with separa¬ 
tions in the range (600,1000) kpc selected to satisfy the 
following constraints: (i) relative radial velocity in the range 
(—175, —75) km s“^ (dark blue histograms; see box in the 
left panel of Fig. [TJ; (ii) tangential velocity in the range 
(0,50) km s“^ (light-blue histogram; see box in the right 
panel of Fig. [T]); and (iii) Hubble velocity fits in the ob¬ 
served range (lightest blue histograms; see box in Fig. 21). 
The histograms in red correspond to the few systems that 
satisfy all three criteria at once. (The scale of the red his¬ 
tograms has been increased in order to make them visible; 
see legends.) 

In these three panels the isolation criteria for select¬ 
ing pairs tightens from top to bottom, and, as a conse¬ 
quence, the total number of selected pairs decreases. The re¬ 
sults, however, are qualitatively similar for all three criteria. 
As discussed before, the radial velocity constraint favours 
high-mass pairs, whereas the tangential velocity constraint 
favours low-mass ones. The Hubble velocity constraint gives 
results intermediate between the other two. 

Very few pairs satisfy all three constraints simultane¬ 
ously. For the case of medium isolation (second panel from 
the top in Fig. Hi only 14 pairs satisfy the three criteria 
simultaneously in the (137 Mpc)^ volume of the MS-II sim¬ 
ulation. The bottom panel of Fig. |4] shows the mass dis¬ 
tribution of pairs that satisfy all three criteria in the case 
of “medium isolation”, the closest to our observed LG con¬ 
figuration (empty red histogram). The mass distribution of 
pairs that satisfy the three criteria simultaneously peaks at 
~ 2 X 10^^ Mq, a value significantly lower than the mass 
suggested by the timing argument. The distribution is quite 
broad, with an rms of 0.4 dex and a full range of values 
that stretches more than a decade; from 2.3 x 10^^ Mq to 
6.1 X 10^^ Mq in the case of medium isolation. 

This finding is one of the main reasons guiding our 
choice of mass for the halo pairs that we select for resim¬ 
ulation in the APOSTLE project; their distribution is shown 
by the solid black histogram in the bottom panel of Fig. 2) 
We describe the selection procedure of these 12 pairs in de¬ 
tail next. 


3 THE APOSTLE SIMULATIONS 

The APOSTLE project consists of a suite of high- 
resolution cosmological hydrodynamical simulations of 12 
Local Group-like environments selected from large cosmo¬ 
logical volumes of a ACDM universe. Each of these volumes 
has been simulated at three different resolutions with the 
same code as was used for the EAG LE simulation project 
llSchave et al.ll2015l : iGrain et al.ll2015|j . 

Results from a subse t of these simulations hav e al¬ 
ready been presented in ISawala et al.l ||2014| . l2015h and 
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Figure 4. Total mass distributions of all pairs in the Millennium- 
II Simulation (MS-II) with separation and relative radial velocity 
similar to the MW-M31 pair, selected assuming different isola¬ 
tion criteria. The dark blue histogram shows pairs that satisfy 
the distance-radial velocity criterion shown by the box in the left 
panel of Fig.^ The lighter blue histogram shows those pairs satis¬ 
fying the constraint shown by the box in the distance-tangential 
velocity panel of Fig. ^ The lightest blue histogram identifies 
pairs satisfying the “Hubble flow” criterion shown by the box in 
Fig.E] The red histogram corresponds to the intersection of all 
three criteria; their scale has been increased to make them more 
easily visible. From top to bottom, the first three panels show 
results for different levels of isolation, as labelled. The numbers 
in brackets indicate the number of pairs in each histogram. Ar¬ 
rows indicate the median of each distribution. The bottom panel 
shows in red the results for medium isolation (empty histogram), 
and compares them with the mass distribution of the 12 candi¬ 
date pairs selected from the DOVE simulation for the APOSTLE 
simulation project (in solid black). 



Figure 5. Relative radial (squares) and tangential (circles) veloc¬ 
ities vs separation for all 12 candidates selected for resimulation 
in the APOSTLE project. The corresponding values for the MW- 
M31 pair are shown by the starred symbols. The shaded area 
indicate the region where we would expect the radial velocities of 
all 12 candidates to lie if the predictions of the timing argument 
held. 


I Oman et al.l (|2015l 'l. Here we describe the selection proce¬ 
dure for the resimulation volumes together with the physical 
processes included in the simulations. We assess the viabil¬ 
ity of these candidates by analyzing, in Sec.jd] their satellite 
populations and comparing them with observations of the 
satellites of the MW and M31 galaxies. 


3.1 The code 

The APOSTLE simulations were run using a high ly modi¬ 
fied v ersion of the N-body/SPH code, P-Gadget3 (ISoringell 
l2005h . The code includes sub-grid prescriptions for star for¬ 
mation, feedback from evolving stars and supernovae, metal 
enrichment, cosmic reionization, and the formation and en¬ 
ergy output from supermassive black holes/active galac¬ 
tic nuclei (AGN). The details of the code and the sub¬ 
grid p hysics model are d escrib ed in detail in ISchave et al.l 
feoish and iGrain et al.l (l2015tl . When used to evolve a 
cosmologically significant volume, this code gives a galaxy 
stellar mass function in good agreement with observations 
in the rang e of galaxy stellar masses between ~ 10® and 
- 10^^ Mr^^ dSchave et all boish . The APOSTLE runs use 
the same pa r amete r choices as the “reference” model in 
ISchave et al.l (120151') . As discussed below, we find good nu¬ 
merical convergence in the galaxy properties we investigate, 
without further calibration. 


3.2 Candidate selection 

The pairs selected for resimulation in the APOSTLE project 
were drawn from the DOVifl cosmological N-body simula- 

® DOVE simulation is also known as COLOR. 
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Table 2. The parameters of the APOSTLE resimulations. The first two columns list labels identifying each run. The following columns list 
the virial masses of each of the primaries at 2 = 0; their relative separation, radial velocity, and tangential velocity in the DOVE simulation, 
as well as the initial baryonic mass per particle in the hydrodynamical runs. The dark matter particle mass is = (l//bar ~ 1) ^gas, 

where /baj. is the universal baryon fraction. (Dark matter-only runs have a particle mass equal to the sum of rugas + '^dmO The last 
column lists the value of the Plummer-equivalent gravitational softening, which is comoving at early times, but fixed at the listed value 
after 2 = 3. 


Name 

Run 

(resolution) 

■^200 

[1012(Mo)] 

[1012 (Mq)] 

separation 

[kpc] 

Vr 

[ km s“^] 

Ft 

[ km s“^] 

TTlgas 

[101(Mq)] 

^max 

[pc] 

AP-l 

HR/MR/LR 

1.66 

1.10 

850 

-51 

35 

0.99/12.0/147 

134/307/711 

AP-2 

MR/LR 

0.85 

0.83 

809 

-39 

97 

12.5/147 

307/711 

AP-3 

MR/LR 

1.52 

1.22 

920 

-35 

84 

12.5/147 

307/711 

AP-4 

HR/MR/LR 

1.38 

1.35 

790 

-59 

24 

0.49/12.2/147 

134/307/711 

AP-5 

MR/LR 

0.93 

0.87 

828 

-33 

101 

12.5/147 

307/711 

AP-6 

MR/LR 

2.36 

1.21 

950 

-18 

60 

12.7/137 

307/711 

AP-7 

MR/LR 

1.88 

1.09 

664 

-174 

24 

11.3/134 

307/711 

AP-8 

MR/LR 

1.72 

0.65 

817 

-120 

96 

11.0/137 

307/711 

AP-9 

MR/LR 

0.96 

0.68 

814 

-28 

48 

10.9/138 

307/711 

AP-10 

MR/LR 

1.46 

0.87 

721 

-63 

48 

11.0/146 

307/711 

AP-11 

MR/LR 

0.99 

0.80 

770 

-124 

22 

11.1/153 

307/711 

AP-12 

MR/LR 

1.11 

0.58 

635 

-53 

50 

10.9/138 

307/711 


tion described bv ljenkini (j2Q13[ ). DOVE evolved a periodic 
box 100 Mpc on a side assuming cosmological parameters 
consistent with the WMAP-7 estimates and summarized in 
Table m DOVE has 134 times better mass resolution than 
MS-I, with a particle mass of 8.8 x 10® Mq (comparable to 
MS-II). Halo pairs were identified in DOVE using the proce¬ 
dure described in Sec. Ea 

Guided by the discussion in Sec. [2] we chose 12 different 
pairs from the Medlso sample that satisfied, at z = 0, the 
following conditions: 

• separation between 600 and 1000 kpc, 

• relative radial velocity, 14, in the range 
(-250,0) km s-\ 

• relative tangential velocity, 14, less than 100 km s“^ , 

• recession velocities of outer LG members in the range 
defined by the box in Fig. 0 

• total pair mass (i.e., the sum of the virial masses 
of the two primary halos) in the range log(Mtot/MQ) = 
[ 12 . 2 , 12 . 6 ]. 

The relative velocities and separations of the 12 pairs 
are shown in Fig. [5] where each pair is labelled with a num¬ 
ber from 1 to 12 for future reference. The main properties 
of these pairs are listed in Table E and the histogram of 
their total masses is shown (in black) in the bottom panel 
of Fig. [4] The parameters of fits to the recession velocities of 
the onter members are shown by the open starred symbols 
in Fig. O 

Note that the 12 pairs chosen span a relatively small 
range of masses compared to what is allowed according to 
the kinematic constraints described in the previous section. 
The lowest mass of the pairs is 1.6 x 10^^ Mq and the most 
massive pair is 3.6 x 10^^ Mq, with an average mass of 2.3 x 
10^^ Mq and an rms of only 30%. The small mass range of 
the pairs was chosen in order to explore the cosmic variance 
at fixed mass. Given the large range of allowed masses it 
would have been impossible to cover the whole range with 
high resolution resimulations. 

It is interesting to compare the masses of the pairs with 


those estimated from the timing argument applied to each 
pair. This is shown by the shaded region of Fig. [5] which 
brackets the region in distance-radial velocity plane that 
should be occupied by the 12 candidates we selected, given 
their total masses. The timing argument does reasonably 
well in 8 out of the 12 cases, but four pairs fall outside the 
predicted region. Three pairs, in particular, have approach 
velocities as high or higher than observed for the MW-M31 
pair (pairs 7, 8, and 11), demonstrating that our mass choice 
does not preclude relatively high approach velocities in some 
cases. 

3.3 Resimulation runs 

Volumes around the 12 selected pairs were carefully config¬ 
ured in the initial conditions of the DOVE simulation so that 
the inner 2-3 Mpc spherical regions centred on the mid-point 
of the pairs contain no “boundary” part icles at z = 0 (for 
details of this “zoom-in” technique, see [Power et al.l l2003l : 
Ijenkind l2013l l. For each candidate pair, initial conditions 
were constructed at three levels of res olution using s econd- 
order Lagrangian perturbation theory ilJenkinsll201Cf) . (Fur¬ 
ther details are provided in the Appendix A.) Resolution 
levels are labelled LR, MR, and HR (low, medium, and 
high resolution), where LR is equivalent to DOVE, MR is 
a factor of ~ 10 improvement in mass resolution over LR, 
and HR is another factor of ~ 10 improvement over MR. 
The method use d to make the zoom initial conditions has 
been described in ijenkind (l2013l ). Particles masses, softening 
lengths, and other numerical parameters for each resolution 
level are summarized in Table E 

Each run was performed twice, with one run neglecting 
the baryonic component (i.e., assigning all the matter to 
the dark matter component, hereafter referred to as dark- 
matter-only or “DM0” runs) and the other using the EAGLE 
code described in Sec. Note that the cosmology assumed 
for the APOSTLE runs is the same as that of DOVE (see 
Table [T}, which differs slightly from that of EAGLE. These 
differences, however, are very small, and are expected to 
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have a negligible effect on our results. At the time of writing, 
only two volumes had been completed at the highest (HR) 
resolution; AP-1 and AP-4. 


4 LOCAL GROUP SATELLITES 

In this section we explore the properties of the satellite pop¬ 
ulation that surrounds each of the two primary galaxies of 
the APOSTLE resimulations to assess whether our choice of 
total mass and subgrid physics gives results consistent with 
the broad properties of the satellite systems of the MW and 
M31. We examine consistency with properties that may be 
considered problematic, given that our choice for the to¬ 
tal mass is significantly lower than indicated by the timing 
argument. In particular, we focus on (i) whether the total 
number of satellites brighter than Mv ~ —8 and Mv ~ —9.5 
(or, equivalently, stellar masfl, Mgai, greater than 1.4 x 10® 
or 5.6 X 10® Mq, respectively), is comparable to the observed 
numbers around the MW and M31; (ii) whether satellites as 
massive/bright as the Large Magellanic Cloud or M33 are 
present; (iii) whether the velocity dispersion of the satellites 
is consistent with observations; and, finally, (iv) whether the 
observed satellites within 300 kpc of either primary (i.e., the 
MW or M31), are gravitationally bound given our choice of 
total mass for the LG candidates. 


4.1 Satellite masses/luminosities 

4-1.1 Number of bright satellites 


The total number of satellites above a certain mass is ex¬ 
pected to be a sensitive func tion of the virial masses of 
the primary hal os (see, e.g., iBovlan-Kolchin et al.l l201ll : 
IWang et al.ll20l3 ~). so we begin by considering the number of 
satellites whose stellar masses exceed 1.4 x 10® M© or, equiv- 
alentl}0, that are brighter than Mv = —8, as a function of 
the virial mass ,M 2 oo, of each host. This brightness limit 
corresponds roughly to the faintest of the “classical” dwarf 
spheroidal companions of the MW, such as Draco or Ursa 
Minor. We choose this limit because there is widespread con¬ 
sensus that surveys of the surroundings of the MW and M31 
are complete down to that limit (see, e.g.. IWhiting et al.l 
I2007I : iMcConnachie et ^l2009l : iMcConnachig 2012l 'l. 

We compare in the top panel of Fig. [B] the number of 
simulated satellites within 300 kpc of each primary galaxy 
to the numbers observed around the MW and M31 (shown 
as horizontal line segments). Satellite numbers correlate 
strongly with virial mass, as expected, and it is encourag¬ 
ing that the observed number of satellites of the MW and 
M31 (shown by short horizontal line segments in Fig. are 
well within the range spanned by our simulations. Since a 


^ Stellar masses are measured within a radius, equal to 15% 
of the virial radius, r200i of the surrounding halo. For satellites 
we estimate rg^i from its peak circular velocity, Umax, and the 
relation between Umax and r2oo for isolated galaxies. 

® We assume a constant mass-to-light ratio of 1 in solar units 
for converting stellar masses to U-band luminosities. This is done 
only for simplicity, since our primary aim is to assess consistency 
rather than to provide quantitative predictions. Later work will 
us e spectrophotometric models and internal extinction as laid out 
in iTravford et all 1I2OI5II . 



log(M2oo/M3,,) 

Figure 6. Top: Number of satellites with stellar mass greater 
than 1.4 X 10^ Mq (or brighter than My = —8) found within 300 
kpc of each of the primary galaxies in our APOSTLE runs plot¬ 
ted versus the primary’s virial mass. We show results for medium 
resolution (MR, blue plus symbols) and high resolution (HR, red 
crosses). The thin dashed lines connects the HR halos to their 
counterparts in MR. The solid, coloured lines indicate the best 
fit to MR data with unit slope; the shaded areas indicate the rms 
range about this fit. Small horizontal lines indicate the observed 
numbers for the MW and M31. Arrows indicate virial masses 
estimated for th e MW and M31 from abundance matching of 
IGuo et~^ ]l l2010ll . Th e “error bar” around each arrow spans dif¬ 
ferent predictions by iBehroozi et al.l ll2013l l and iKravtsov et al.l 
<l2014ll . Bottom-. Same as the top panel but for satellites with stel¬ 
lar mass greater than 5.6x 10® Mq (or brighter than Mv = —9.5). 


lower-mass limit of Mgai = 1.4 x 10® Mq corresponds to just 
a few particles at MR resolution, we repeat the exercise in 
the bottom panel of Fig. (6] but for a higher mass limit; i.e., 
A/gai > 5.6 X 10® Mq {Mv < —9.5). The results in either 
panel are reassuringly similar. 

The impact of numerical resolution may be seen by com¬ 
paring the results for the medium resolution (MR; blue “-I-” 
symbols) with those obtained for the HR (high-resolution) 
runs in Fig. |6] (red crosses). On average, the number of satel¬ 
lites increases by only ~ 10% when increasing the mass reso- 
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Figure 7. Stellar mass of the primary galaxies (circles) and 
their brightest satellites (crosses) for our 12 APOSTLE simula¬ 
tions, as a function of their maximum circular velocity. The large 
black solid symbols denote observational estimates for the two 
brightest members of the Local Group and their brightest satel¬ 
lites: MW, M31, LMC, and M33, respectively. Grey curves ind i¬ 
cate “abundan c e ma t ching” predictions fro m IGuo et alj ll201Cl) : 
iBehroozi et alJ ll2ni3h : rKravtsov et al.l<l2014h . Note that satellites 
as bright and massive as the LMC and M33 are not uncommon 
in our simulations. Note also that our primaries tend to be un- 
dermassive relative to observations and to predictions from abun¬ 
dance matching. The reverse applies at lower circular velocities. 
Black dashed lines connect matching systems at different resolu¬ 
tion levels. See text for further discussion. 


lution by a factor of 10, indicating rather good convergence. 
One of the halos in HR, however, hosts almost twice as many 
satellites as its counterpart in MR. The reason is that a rel¬ 
atively large group of satellites has just crossed inside the 
300 kpc boundary of the halo in HR, but is still outside 300 
kpc in MR. (We do not consider the low-resolution LR runs 
in this plot because at that resolution the mass per particle 
is 1.5 X 10® Mq and satellites fainter than My ~ —13 are 
not resolved, see Table [2]) 

We also note that, had we chosen larger masses for our 
LG primaries, our simulations would have likely formed a 
much larger number of bright satellites than observed. For 
example, for a virial mass of ~ 2 x lO'^^ Mq (the value esti¬ 
mated for the MW by abundance-matching analyses; see the 
upward arrow in Fig. El) our simulations give, on average, 25 
satellites brighter than My = —8 and 15 satellites brighter 
than My = —9.5. These are well in excess of the 12 and 
8 satellites, respectively, found in the halo of the MW. The 
same conclusion applies to M31, where abundance matching 
suggests a halo mass of order 7 x 10^^ Mq. It is clear from 
Fig. El that our simulations would have produced a number 
of satellites well in excess of that observed around M31 for 
a halo as massive as that. 

Satellite numbers are, of course, sensitive to our choice 
of galaxy formation model, and it would be possible, in prin¬ 
ciple, to reduce the number of bright satellites by reducing 
the overall galaxy formation efficiency in low-mass halos. It 



1^200 [l^sun] 

Figure 8. Top-. Radial velocity dispersion of luminous satellites 
within 300 kpc of each primary in the APOSTLE runs, as a func¬ 
tion of the host’s virial mass. The red crosses and blue plus signs 
correspond to the resolution levels, HR and MR, respectively. 
The thin dashed lines connects halos of different resolutions. The 
solid coloured line indicates the best power-law fit cr oc the 

shaded areas indicate the rms around the fit. Note the strong cor¬ 
relation between velocity dispersion and host mass; the relatively 
low radial velocity dispersion of the MW satellites is best accom¬ 
modated with a fairly low virial mass, lower than expected from 
abundance matching (see upward pointing arrow, the same as in 
Fig. EJ. Bottom: Same as top, but for the line-of-sight velocity 
dispersion of one satellite system as seen from the other primary. 
Compared with the observed line-of-sight velocity dispersion of 
M31 satellites, shown by the horizontal line, this measure favours 
a lower virial mass than inferred from abundance-matching anal¬ 
yses, indicated by the arrow. 


is interesting that the same galaxy formation model able to 
reproduce the g alaxy stellar mass function in large cosmo¬ 
logical volumes (|Schave et al.|[2015l) can reproduce, without 
further tuning, the number of satellites in the APOSTLE 
resimulations. 
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Figure 9. Top: Radial velocity vs distance for all satellites in 
the medium-resolution (MR) APOSTLE runs (“plus” symbols), 
compared with the observed radial velocities of the MW luminous 
satellites (filled squares). The positions of observed satellites in 
this phase-space correspond to regions that are well populated 
in the simulations. Very few satellites are “unbound” (marked 
by circles), as judged by escape velocity curves computed assum¬ 
ing a Navarro-Frenk-White profile of given virial mass (see dot¬ 
ted lines). Bottom: Same as top but for line-of-sight velocities of 
satellites as seen from the other primary. This may be directly 
compared with data for M31 satellites (filled triangles). 


4.1.2 Most massive satellites 

Another concern when adopting a relatively low virial mass 
for the LG primaries, is that satellites as bright and mas¬ 
sive as M33 in the case of M31, or the LMC in the case of 
the MW, may fail to form. We explore this in Fig. [7] where 
we plot the stellar mass of the primaries (open circles) and 
that of their most massive satellite found within 300 kpc 
(crosses, “plus”, and square symbols for HR, MR, and LR, 
respectively) as a function of their peak circular velocity. 
We also plot, for reference, the rotation velocity and stel¬ 
lar masses of the MW and M31, as well as those of their 
brightest satellites, the LMC and M33, respectively. For the 
LMC, the error bar indicates the velocity range spanned by 


two different estimates, 64 km s ^ from 


ll2002l) . and 87 km s ^ from lOlsen et al 


van der Marel et al.l 


(120111) . For M33, we 


use a rotation speed of 110 km s“^, with an error bar that 
indicat es the range of veloci ties observed between 5 and 15 
kpc bv ICorbelli et al.l ll2014h : we use this range as an esti¬ 
mate of the uncertainty because the gaseous disk of M33 
has a very strong warp in the outer regions which hinders a 
proper determination of its asymptotic circular speed. 

Fig. 0 makes clear that there is no shortage of massive 
satellites in the APOSTLE simulations: 6 out of 24 primaries 
in the MR runs have a massive subhalo with Vmax exceeding 
100 km s“^ (comparable to M33), and 12 out of 24 have 
Knax > 60 km s“^ (comparable to the LMC). This result is 
robust to numerical resolution; the numbers above change 
only to 7 and 11, respectively, when considering the lower- 
resolution LR runs (see open squares in Fig. [7]|. 

With hindsight, th i s resu lt is not entirely surprising. Ac¬ 
cording to I Wang et al.l (l2012l) . the average number of satel¬ 
lites scales as 10.2 (i//0.15)“^'^^, where v = Fmax/Vioo is 
the ratio between the maximum circular velocity of a sub¬ 
halo and the virial velocity of the primary halo. A halo of 
virial mass 10^^ Mq has V 2 Q 0 = 145 km s“^ , so we expect, 
on average, that 1 in 11 of such halos should host a satel¬ 
lite as massive as M33 and 1 in 2.3 one like the LMC0. Our 
results seem quite consistent with this expectation. 

A second point to note from Fig. [7] is that the stellar 
mass of the primary galaxies are below recent estimates for 


eithe r the MW (5.2 x 10^“ M©, according tojBov^^^ix 


2013) or M31 (1.03 x 10^^ Mq, according to iHammer et al 


20071) . Indeed, none of our 24 primaries have masses exceed¬ 
ing 5 X 10^° Mq, a consequence of the low virial masses of our 
selected LG pairs coupled with the relatively low galaxy for¬ 
mation efficiency of the EAGLE “Ref” model in ~ 10^^ Mq 
halos. This may be seen in Fig. 0 where the open circles he 
systematically below the grey lines, which show abundance 
matching predictions taken from three recent paper s. _ 

This issue has been discussed bv ISchave et al.l ll2015h . 
and is reflected in their fig. 4, which shows that the EA¬ 
GLE “Ref” model we use here underpredicts the number of 
galaxies with stellar masses a few times 10^° Mq . That same 
figure shows that the opposite is true for smaller galaxies: 
the “Ref” model actually overpredicts slightly but systemat¬ 
ically the number of galaxies with stellar masses a few times 
10® Mq. 

Galaxy formation efficiency in our runs thus seems 
slightly too low in MW-sized halos, and slightly too high 
in systems of much lower mass. This slight mismatch in the 
galaxy mass-halo mass relation manifests itself more clearly 
on LG scales, making it difficult to match simultaneously the 
stellar masses of the MW and M31 as well as that of their 
luminous satellites. Indeed, increasing the halo masses of our 
LG candidates would lead to a better match to the stellar 
masses of the primaries, but at the expense of overpredicting 
the number of bright satellites (see, e.g.. Fig. O, unless the 
star formation and feedback parameters are recalibrated. 

The stellar masses of LG galaxies are therefore a sensi¬ 
tive probe of the galaxy formation efficiency on cosmological 
scales, and provide an important constraint on the ability of 


° These numbers assume a maximum circular velocity of 
100 km for M33 and 60 km s“^ for the LMG. 
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cosmological codes to reproduce the observed galaxy popu¬ 
lation. 

4.2 Satellite kinematics 

We can also use the kinematics of MW and M31 satellites to 
gauge the consistency of our results with the simulated satel¬ 
lite population. We begin by considering, in the top panel 
of Fig. |8] the radial velocity dispersion, Or, of all satellites 
more massive than 1.4 x 10 ® Mq (brighter than Mv = —8) 
within 300 kpc from either primary. We expect the satellite 
velocity dispersion to scale as ar oc V 200 oc Mjgg, so the 
solid line indicates the best fit with that slope to the data 
for our 24 systems. (Symbols are as in Fig. [b] blue -|- sym¬ 
bols indicate MR resolution, red crosses correspond to HR 
resolution.) This scaling describes the correlation shown in 
Fig. [8] fairly well, and its best fit suggests that a 10^^ Mq 
system should host a satellite system with a radial velocity 
dispersion of 98 ± 17 km s“^. (The shaded area in Fig. [8] 
indicates the rms scatter about the best fit.) The ob¬ 

served dispersion of the MW satellites (86 km s“^; shown 
by the dashed horizontal line) thus suggests a mass in the 
range 3.9 x 10^^ M© to 1.1 x 10^^ M©. The scatter, however, 
is large and these data alone can hardly be used to rule out 
larger virial masses. 

The bottom panel of Fig. [8] shows a similar analysis, but 
applied to the satellites of M31. Since M31-centric radial ve¬ 
locities are not available, we repeat the analysis for our 24 
satellite systems using projected velocities, measured along 
the line of sight from the other primary of each pair. (The 
analysis uses only satellites within 300 kpc from the galaxy’s 
center.) The relation between the projected velocity disper¬ 
sion, (Jvios, and virial mass is again reasonably well described 
by the expected Mjgg scaling. In this case, the observed pro¬ 
jected dispersion of 99 km s“^ for the M31 satellites suggests 
a mass in the range 7.6 x 10^^ M© to 2.1 x 10^^ Mq although 
the scatter is again large. 

The main conclusion to draw from Fig.[8]is that the ob¬ 
served velocity dispersions of both the MW and M31 satel¬ 
lites are well within the ranges found in our simulations. 
The relatively low mass of our pairs thus does not seem to 
pose any problems reproducing the kinematics of the LG 
satellite population. On the contrary, as was the case for 
satellite numbers, the satellite kinematics would be difficult 
to reconcile with much larger virial masses. For example, if 
the MW had a virial m ass of 1.8 x 10^^ Mq , as suggested 
by abundance-matching iIGuo et al.ll2010l f. then its satellite 
systems would, on average, have a velocity dispersion of or¬ 
der 120 km s“^, which would exceed the 86 km s“^ observed 
for the MW companions. 

Finally, we check whether any of the MW or M31 satel¬ 
lites would be unbound given the mass of the primaries cho¬ 
sen for our sample. We explore this in Fig. [9l where we plot 
the observed velocity (radial in the case of the MW; line-of- 
sight in the case of M31) of satellites as a function of their 
distance to the primary’s center (solid symbols). 

The dotted lines delineate the escape velocity as a func¬ 
tion of distance for halos with virial massses of 6 x 10^^ 
and 2 x 10^^ Mq , corresponding to roughly the minimum 
and maximum virial masses of all primaries in our sample. 
The escape velocities assum e a Navarro-Frenk-White profile 
dNavarro et al.l [19961 . ll997l L with a concentration c = 10. 


“Escapers” (i.e., satellites with 3D velocities exceeding the 
escape speed for its primary) are shown by circled symbols. 
These are rare; only 3 of the 439 Mv < —8 satellites exam¬ 
ined are moving with velocities exceed ing the nominal NFW 
escap e velocity of their halo (see also iBovlan-Kolchin et al.l 

l20ial . 

It is clear from this figure that, for our choice of masses, 
none of the MW or M31 satellites would be unbound given 
their radial velocity. Indeed, even Leo I and And XIV (the 
least bound satellites of the MW and M31, respectively) are 
both within the bound region, and, furthermore, in a region 
of phase space shared with many satellites in our APOSTLE 
sample. The kinematics of the satellite populations of MW 
and M31 thus seems consistent with that of our simulated 
satellite populations. 


5 SUMMARY AND CONCLUSIONS 

We have analyzed the constraints placed on the mass of the 
Local Group by the kinematics of the MW-M31 pair and 
of other LG members. We used these constraints to guide 
the selection, from a large cosmological simulation, of 12 
candidate environments for the EAGLE-APOSTLE project, 
a suite of hydrodynamical resimulations run at various nu¬ 
merical resolution levels (reaching ~ 10^ Mq per gas particle 
at the highest level) and aimed at studying the formation of 
galaxies in the local Universe. 

APOSTLE uses the same code and star forma¬ 
tion/feedback subgrid modules developed for the EAGLE 
project, which yield, in cosmologically-representative vol¬ 
umes, a galaxy stellar mass function and average galaxy sizes 
in good agreement with observations. This ensures that any 
success of our simulations in reproducing Local Group-scale 
observations does not come at the expense of subgrid module 
choices that might fail to reproduce the galaxy population at 
large. We also compare the simulated satellite populations of 
the two main galaxies in the APOSTLE resimulations with 
the observed satellite systems of M31 and MW to assess 
consistency with observation. 

Our main conclusions may be summarized as follows: 

• The kinematics of the MW-M31 pair and of other LG 
members are consistent with a wide range of virial masses for 
the MW and M31. Compared with halo pairs selected from 
the Millennium Simulations, the relatively fast approach ve¬ 
locity of MW and M31 favours a fairly large total mass, of 
order 5 x 10^^ Mq. On the other hand, the small tangential 
velocity and the small deceleration from the Hubble flow of 
outer LG members argue for a significantly smaller mass, of 
order 6 x 10^^ Mq. Systems that satisfy the three criteria are 
rare—only 14 are found in a (137 Mpc)® volume—and span 
a wide range of masses, from 2.3 x 10^^ to 6.1 x 10^^ Mq, 
with a median mass ~ 2 x 10^^ Mq . 

• Given the wide range of total masses allowed, the 12 
candidate pairs selected for resimulation in the APOSTLE 
project were chosen to loosely match the LG kinematic 
criteria and to span a relatively narrow range of masses 
(from 1.6 X 10^^ to 3.6 x 10^^ Mq, with a median mass of 
2.3 X 10^^ Mq). This enables us to explore the cosmic vari¬ 
ance of our results at fixed mass, and, potentially, to scale 
them to other mass choices, if needed. 
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• Large satellites such as LMC and M33 are fairly com¬ 
mon around our simulated galaxies, although their total 
virial mass is well below that estimated from the timing 
argument. 

• The overall abundance of simulated satellites brighter 
than Mv = — 8 is a strong function of the virial mass as¬ 
sumed for the LG primary galaxies in our simulations. The 
relatively few (12) such satellites around the MW suggests a 
fairly low mass (~ 6x 10^^ Mq); the same argument suggests 
a mass for M31 about twice as large (~ 1.2 x 10^^ Mq). 

• The velocity dispersions of simulated satellites are con¬ 
sistent with those of the MW and M31. This diagnostic also 
suggests that virial masses much larger than those adopted 
for the APOSTLE project would be difficult to reconcile with 
the relatively low radial velocity dispersion observed for the 
MW satellite population, as well as with the projected ve¬ 
locity dispersion of the M31 satellite population. 

• The primary galaxies in the simulations are less mas¬ 
sive than current estimates for the MW and M31. The most 
likely reason for this is an inaccuracy in the subgrid mod¬ 
eling of star formation and feedback and its dependence on 
halo mass. 

Our overall conclusion is that, despite some shortcom¬ 
ings, the APOSTLE simulation suite should prove a wonder¬ 
ful tool to study the formation of the galaxies that populate 
our cosmic backyard. 
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APPENDIX A: PARAMETERS OF THE 
INITIAL CONDITIONS 


The initial conditions for the DOVE and APOSTLE simula- 
tions were gen erated from the PANPHASIA white-noise field 
lljenkind 20131) using second-order Lagrangian perturbation 
theory jjenkinsll201(]l l. The coordinates of the centers and 
the radii of the high resolution Lagrangian regions in the 
initial conditions, as well as the positions of the MW and 
M31 analogs at 2 = 0, for the twelve APOSTLE volumes are 
given in Table \M\ 
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Table Al. The positions of main halos at 2 : = 0 and parameters of the high resolution Lagrangian regions of the APOSTLE volumes 
in the initial conditions. The first column labels each volume. The next columns list the (X,Y,Z) coordinates of each of the primaries at 
2 : = 0, followed by positions of the centers and the radii of the Lagrangian spheres. The phase descriptor for the APOSTLE runs is, in 
PANPHASIA format, [Panphl,L16,(31250,23438,39063),S12,CH1292987594,DOVE]. 


Name 

[Mpc] 

Yi 

[Mpc] 

Zi 

[Mpc] 

^2 

[Mpc] 

Y 2 

[Mpc] 

Z 2 

[Mpc] 

Xi 

[Mpc] 

Yi 

[Mpc] 

Zi 

[Mpc] 

H, 

[Mpc] 

AP-l 

19.326 

40.284 

46.508 

18.917 

39.725 

47.001 

26.5 

39.1 

39.0 

5.7 

AP-2 

28.798 

65.944 

17.153 

28.366 

65.981 

16.470 

28.1 

60.2 

18.4 

10.2 

AP-3 

51.604 

28.999 

11.953 

51.091 

28.243 

12.061 

46.0 

31.7 

11.6 

12.0 

AP-4 

63.668 

19.537 

72.411 

63.158 

20.137 

72.467 

57.1 

20.6 

74.9 

6.0 

AP-5 

42.716 

87.781 

93.252 

42.872 

88.478 

93.671 

40.8 

85.4 

91.8 

9.9 

AP-6 

35.968 

9.980 

43.782 

36.171 

9.223 

43.251 

32.9 

13.1 

45.2 

7.1 

AP-7 

91.590 

43.942 

14.826 

91.822 

43.323 

14.885 

99.3 

39.7 

15.9 

7.1 

AP-8 

4.619 

22.762 

85.535 

4.604 

23.508 

85.203 

4.9 

20.4 

89.9 

9.3 

AP-9 

57.044 

88.490 

74.765 

57.496 

87.889 

74.456 

55.2 

93.4 

76.5 

5.7 

AP-10 

61.949 

24.232 

98.305 

61.867 

24.925 

98.124 

62.5 

24.5 

93.5 

8.4 

AP-11 

12.564 

48.080 

35.249 

12.484 

47.793 

35.959 

18.3 

43.1 

29.9 

5.9 

AP-12 

97.553 

89.587 

72.093 

97.351 

90.100 

72.407 

98.5 

91.9 

81.9 

5.7 





